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I Abstract. Sifting limits for the A'^A" sieve, Selberg's lower bound sieve, are 

computed for integral dimensions 1 < k < 10. The evidence strongly suggests 
that for all K > 3 the A^A~ sieve is superior to the competing combinatorial 
sieves of Diamond, Halberstam, and Richert. A method initiated by Grupp 
and Richert for computing sieve functions for integral k is also outlined. 

Q 

^ 1. Introduction 

^ Let A he a. sequence of integers, and V a set of primes. Recall that the goal of 

the sieve method is to obtain bounds for 

^ (1) SiA,V,z)= Yl 1' 

where 

^ P{z) =11 p. 
, pev 

One expects that the sequence A is well-behaved in that Ad, the elements of A 



> 

divisible by d, satisfy 
00 X 

cn (2) l-^di = -^ + 7^d, 

where f{d) is some multiplicative function, and the errors, TZd, are relatively small, 
at least on average. In fact, suppose that there exists a constant A>1 such that 

>< """^ 

^ In addition, one assumes that 



(4) ^l^ = ,iog. + 0(l), 

p<s ■' 

and refers to k as the dimension, or density, of the sieve. 



H. Diamond and H. Halberstam, in association with the late professor H.-E. 
Richert, constructed a class of sieves for all dimensions k > 1. Their sieves (DHR 
sieves for short) combine elements of Selberg's A^ upper bound sieve and the combi- 
natorial sieves of Rosser-Iwaniec. For an account of their work, we refer the reader 
to their recent book An important parameter in a sieve is the sifting limit (3^^, 
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beyond which the lower bound sieve yields a positive lower bound. The calculations 
in Chapter 17 of [1] show that for the DHR sieves, (3^ < 2.44k. 

Selberg investigated an alternative lower bound sieve method, known as the 
A^A~ sieve, for large dimensions k. The starting point for this sieve, similar to 
the upper bound sieve, is the observation that for any set of real numbers A^, 
normalized so that Ai = 1, 

2 



s{A,r,z) > ^ 








neA 


p\n 

\ P<z / 




d\n 

\d\P{z) 1 



Selberg proved that for sufficiently large k, this sieve yields /3k ^ 2k + 19/36. As a 
consequence, the A^A~ sieve is superior to the DHR sieves if k is taken sufficiently 
large. How large is sufficiently large? For small integer k with 2 < k < 10, we 
prove 

Theorem 1. Suppose S {A,V, z) is as defined in ([!]), and that A satisfies (|3|, 
and Q. Letting \A\ — x, and z — x^^^" , we have 

S{A,V,z)-» 
for pairs k and listed in the table below. 



log'^ X 





2 


3 


4 


5 


6 


7 


8 


9 


10 


(3. 


4.516 


6.520 


8.522 


10.523 


12.524 


14.524 


16.524 


18.525 


20.525 



Thus, Selberg's sifting limit is approached rapidly from below. Indeed, although 
we have restricted the argument to integer 2 < k < 10, we expect that /3„ < 
2k + 19/36 for all k. When compared with the DHR sieves, the A^A~ sieve gives 
a better sifting limit for integral k > 3. The table below gives a comparison of 
the two sieves. 

Table 1. Sifting Limit Comparison 



K 


2 


3 


4 


5 


6 


7 


8 


9 


10 


DHR /?« 


4.266 


6.640 


9.072 


11.534 


14.014 


16.504 


18.998 


21.495 


23.992 


A^A- 


4.516 


6.520 


8.522 


10.523 


12.524 


14.524 


16.524 


18.525 


20.525 



More improvements are certainly possible. Recently, Sara Blight [U pp. 28-29] 
has shown that /32 < 4.45, /^a < 6.458, and (34 < 8.47. Her work features a set 
of weights that take into account numbers composed of up to three prime factors. 
These weights were suggested by Selberg as a modification to the A^A~ sieve. 

One interesting application of these sieves is to almost-primes in polynomial 
sequences. In a forthcoming paper, the author will show that a weighted A^A~ 
sieve is capable of producing better results than the weighted DHR sieves when 
the polynomial is a product of linear irreducible factors, for example. However, 
the DHR sieves still perform quite well in the higher dimensional setting when the 
irreducible factors of the polynomial are each of a large degree, owing to the optimal 
nature of the DHR construction when k = 1. 
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2. Sieve Setup 



Following Selberg, we define f := f * fi and let be an arbitrary sequence of 
real numbers with the property that = if d is not squarefree, or if d > ^. Next, 
define by the relation 



(5) 



fir) 



E 



A, 



dr 



f{dr) 



By Mobius inversion, we also have 

^^{d)\d 
fid) 



y- 



f'idr) 



E i^E 



(n,P(2)) = l 



neA 



are 


constant 


Ad 
Ai 


. Since A 


/ 


\ 


1 


-El 


V 





E 



The right-hand side can be rearranged using a well-known identity. In particular, 
we have 



Lemma 1. With defined as in we have 



(6) 

where 
(7) 

and 
(8) 



E 

neA 



1 ] 




( ] 


2 


E 




E 


= |^|6^ + 


din 

\d\Piz) / 




i/|n 





e.-E E 

m d 

(d,m) = l 



/'(m) f{d) 



Ea^wc-- 



d,lii ,V2\P{z) 



For our purposes, we divide both sides of this identity by \\ and choose 



(9) 



1, ifd=l, 
= ^ — 1, if c? is prime and d < z, 
0, otherwise. 



The identity in ^ distinguishes ([t]) as the main term and ([8| as the error term for 
the sum. This identity is the starting point of the A^A~ method and has appeared 
in various forms in the works of Selberg [HI See Section 7 on p. 82], Bombieri [21 
See Theorem 18 on p. 65], Cojocaru and Murty [31 See Theorem 10.1.1 on p. 178], 
Greaves [SJ See Lemma 1 on p. 286], and others. 
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To produce a positive lower bound for ([T]) we will show that 
\A^ + ^>\A\V{z) (c + o(l)), 
where c is some small positive constant, and 

To begin, suppose that |^| = x, and let z — x^l^ . It is easy to see that 

^log^x. 

Next, choosing 

(10) zi^ = x^-\ 

and recalling ([3|, we have 

^ ra<zi^ d,vi,U2 m<z(^ 

m\P{z) [d,ui,U2]=m m\P{z) 

« E m'M7''('") |7^™| « — 
mT^e log ^ 

Here we have used the fact that ^ is bounded, which will be explained below. 
The (r will be chosen as 

where P(w) is a polynomial that is positive for < w < u. Therefore, 

(11) ^i^Effey^ ,«^P E 77L« E T^T^w^- 

r|P(2:) 

In the case when = 1, the are well-understood. We will refer to this choice of 



Xi, as Xi,. It is known, for example, that 



< 



Ai 



Since 



(12) |Ad| < sup |P(u;)|Ai, 

0<w<u 

and 

^2(r) /logC/r 



it is clear that 



r|P(z) 



/'(r) V log^ / Q<w<u 



sup |P(w 



All - inf \P{w) 
It follows that the sequence 

is bounded. 



SIFTING LIMITS FOR THE A^A" SIEVE 



5 



Finally, since 



\A\^ + ^ = \A\V{z) 



6^y(z) 1 (Ba 



we have 



1^1 



A? 



\A\V{z) 



O 



1 

log X 



We showed in (11) that XiV{z) is bounded, and so our priority is in the analysis of 
6aV{z). 



3. Analysis of the Main Term 



In this section, we will treat the expression 6_4 occurring in the main term of the 
A^A~ lower bound sieve. First, let us recall that with Selberg's choice of weights 
cid in ([9| we have that 



/i^(m) 



m\P{z) 



^ /'(m) ^ f(p) 

m\P(z) 



(Cm Cpm) 1 



upon omitting the condition that {m,p) = 1 in ([t]). We wish to smooth this 
expression using the asymptotic formulas for 



(14) 

and 
(15) 



G(r,z) 



^ /'(to) 

m<r ■1^1 
m\P{z) 



■ / \ogr \ 
y log z J 

V{z) 



where jK,{u) is the continuous solution of the differential delay equation 
(16) uj\u) = Kj{u) - Kj{u - 1), 

for u > 1, with 

e — ''^ 



(17) 



if < u < 1, 



r(K + i) 

0, if u < 0. 



We remark that if k is held fixed, then Jk(u) increases to 1. Now, using Riemann- 
Stieltjes integration and replacing the integrators with their corresponding smooth 



approximations in (14) and (15 1, we expect that 
(18) 

6a > 



1 .2,- /^logr 



V{z) A 



logz/ V{z) 



1 Ji 



(Cr - Qsr) ^ZZ—d-Jn 



logs 



logz , 
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This is indeed the case since, more specificaUy, if one regards k and u :— 
as fixed, then one has 



> 1 



(19) 

and 
(20) 



m\P(z) 



fl'^{m) (logx) 

/'(m) = V{z) 



l + O 



1 



\ogz 



p<s ^ 



making the error in (18) of order at most logz)~^. The formula in (20 1 is 

merely our assumed density hypothesis in (|4|. On the other hand, the bound in 
( 19 1 is a consequence of 

Lemma 2. For any t = > 0, we have 



G(r,z) 



V{z) 



O 



logz 



Lemma [2] is discussed in some detail in Halberstam and RichertO See Section 4 
on p. 197]. Now, let us define 

- ^°S^>1, 



and 
(21) 
where 
(22) 



Cr=P 



logO^ 
log z 



logz 
logz 



if r > 



is a polynomial in the range r < ^. Using these definitions simplify the integrals 



occurring in the analysis of &Ai a-nd making the variable change v — and 

log 2 ' 



t = |2££, in m, we have 



6a> 



1 



V{z[ 



where 



and 



P* {u-vff, iv)dv, 



(p* _ _ p* _ „ _ t)f jj'^ (v) dv. 



Furthermore, after making the change of variable w — u — v, and using (22 1, these 
integrals further simplify to 



Ii=l P{wf i'^{u-w)dw, 



and 



^0 



,2 dt 



{P{w)-P* {w^t)Y -j',{u-w)dw. 
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The latter integral must be further dissected to account for the vanishing of the 
function P* in the range w < t < 1. As a result, the region in the latter integral 
naturally splits into three distinct pieces. After splitting the range of integration 
to account for this, we find that 



(23) 




V{z) 




where 








(24) 


= 


r 2 

Ja 

I {P{w)-P{w^t)fjj'Ju- 




(25) 


^2 = 


w) dw, 


(26) 




pi nW 1, 

/ / iP{nj)-P{w~t)f-j'Au- 
Jo Jo ^ 


w) dw, 








(27) 


Ti = 


[ f P [w f '^j'^ (u - w) dw. 

Jo Jw 





Contrary to initial appearances, the innermost integral in (25) and (26 1 docs not 
have a singularity at i = because the constant term does not appear in the 
difference P (w) — P {w — t). The next step is to employ a device of Grupp and 
Richert to evaluate these integrals. Before moving on, let us remark that if u 
K — 1/3 — d, and k is taken sufficiently large, SelbergjSl See pp. 174-176] has shown 
that if one sets P{w) = w + a, one has 



V{z) V 2 18 

Choosing a so that d is as large as possible with —a^ + — > Oi we see 

that the optimal choice is a = 1/4, which implies that a positive lower bound is 
achieved when d < —7/72. A slightly more complicated argument that involves a 
more sophisticated set of weights will give d < —7/72, and this is enough to show 
that the sifting limit > 2m + 1 = 2k + ||, upon taking d — —7/72. The weights 
that achieve this involve divisors of n consisting of two and three prime factors. As 
the author's investigations of the use of higher degree polynomials in this problem 
has not met with much success, we will follow Selberg and restrict our attention to 
linear polynomials as well. 



4. The Kn{u,X) Functions 

In order to evaluate the integrals arising in our sieve, we will need to decompose 
j'l^. In his dissertation, Wheeler See Proposition 3.1.1 on p. 18] noted that j^, 
as well as its derivatives, could be decomposed into a sum of simpler functions 
Kn{u, A), each defined for A > —1 and n > 0. More specifically, we have 

(28) j,iu)^-f^ i-^rKn{u,K). 

The sequence of functions Kn{u, A) is defined by the equations 

(29) Ko{u,X)^u^ u>0. 
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and 
(30) 



u> n> 1. 



We also specify that these functions vanish if u < n, and thus 
(31) i^„(u,A) = u<n. 



To justify that the decomposition in ( 28 ) is vahd, one can verify that the expression 
on the right-hand side of (28) satisfies the delay-differential equation in (16). This 
foUows from the observation that for n > 1, 

uK'^{u, k) — Kn-l{u — 1, k) + KKn{u, k) . 



Upon separating the first term from the decomposition in ( 28 1 , we have 



0<n<u 



l<n<ti 



^WTYf'^''^''^^T(^)^^ (-.)".<(«,«:) 



(k + 1) 



l<n<ii 



i (k + Ij i (k + 1) ^ ^ — ' 



^ ' l<n<u 



which is 



0<n<M 



-K7 

X'^Ti) ^ 



r(K + i) „f^/ ' ^ ' r(K + i) , f-^, ^ ^ " ' 



<n<u 



0<n<u 



(ac + 1) 



0<n<n-l 



Therefore, the expression occurring on the right-hand side of (28) satisfies the same 
differential delay equation as the sieve function Jk(u). We will only be concerned 
with integral dimensions k throughout this discussion, and thus will focus on inte- 
gral A > — 1. In fact, the most important case occurs when A = and the following 
lemma will provide us with a useful tool to understand the cases when A 7^ 0. 

Lemma 3. // A > 1 and n > 0, then 

d 



du 



X„(u,A) = AX„(u,A-l) 



Proof. We assume that u > n, for the result is obvious otherwise. Our proof is by 

3 see that 



induction on n; the case n = is obvious. From (30), we see that 
d 



(32) ^K^+i{u, A) = ^'^^ + Xu'-' 

du u 



n+l 



-dt. 
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On the other hand, we can use (30 1 together with the inductive hypothesis and 
integration by parts to get 

XKn+i[u, X - 1) ^ u / 

Jn+1 ^ 

U Jn+1 t^+^ 



The desired result follows by comparing this with ( 32 ) . □ 
As an application of this lemma, it is easy to deduce that 

(33) = f?^ E i-^rKn{u, K - 1), 

^ ' 0<n<u 

and indeed expressions for higher derivatives of can be obtained with more 
applications of Lemma [3j if desired. 

5. The Case A = 

Grupp and Richert |6j made a close study of if„(u,0), obtaining useful power 
series representations for these functions. Their notation differs from Wheeler's, 
but their results can be translated easily since 

Kn{u,0) = {u+l)I„ + l{u+l). 

In this section and the following one, we shall write Kn{u) in place of Kn{u, 0). We 
can obtain an analytic continuation of the function Kniu) if we define Kn{z) by 
the equations 

Ko{z) = 1 ^z> -1, 

and 

f-Z J, 

Kn{z)^ Kn-i{z~l)- ^z>n-l. 

It is easy to see that Kn{z) is an analytic function for > n — 1 and coincides 
with Kn{u) for real values of u > n. Thus, the power series 

oo 

(34) i^„(u) =^5j(n,c)(u- (7i + c))J 

3=0 

is valid for \u — (n + c)\ < 1 + c and u > n, and any c > 0. Moreover, the constant 
coefficients satisfy 

(35) 6o («, c) = Kn (n + c). 

Now, we have the following recursive formula for the rest of the coefficients bj(n, c), 
where j =/= 0. 

Lemma 4. If j > 1, n > 1, and c > 0, then 

b,in, c) = ^-^^ Y.{~l)%{n - 1, c){n + c)'. 

"'^ '1=0 
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Proof. From (30) and (34), we obtain i^^(w) = Kn-i{u — and 

oo ^ oo 

y^j6,-(n, c)(u - (n + c))^"^ = - 6;(n - 1, c)(u - (n + c))'. 
j=a 1=0 

If |m — (n + c)| < 1 + c and n > 1, then |u — (rt + c)| < n + c and 

{u- (n + c))'' {-!)'' 



fe=0 



(n + c) 



fc+i 



Inserting this last equation into the previous one, we find that 



k+l 



J2jbjin,c){u-in + c)y ' =J2Y. fa\ ^^fc+i ^'('^ ~ ^' ^) + ^))^ 

j=0 k=0 1=0 ^ ' 

The desired result follows by equating coefficients of (u — (n + c))-'"^ on both 
sides. □ 

An alternative form of the recursive formula for the sequence bj{n,c) will also 
be useful for induction arguments to follow. 

Lemma 5. If j > I, n > 1, and c> 0, then 

1 



bj{n, c) 



j{n + c) 



l,c) - {j ~ l)&j_i(n,c)}. 



Proof. First, observe that from ( [34| , since bj{n,c) are precisely the coefficients in 
the power series expansion of Kn (u) centered about u = n + c, we have 



(36) 



bj{n,c) 



i=n+c 



Next, from (30), we see that 



1=0 



Upon dividing both sides of this equation by (j — 1)! and evaluating at m = n + c, 
the formula follows from ( 36 ) . 



□ 



Grupp and Richert [B] gave the useful bound 



(37) 



\bi(n,c)\ < —, 



valid for < c < 5 and j > 2. We will need a bound in a larger range of c for our 
purposes. Also, we will be content to accept a slightly worse bound in exchange for 
a simpler proof. Thus, we prove 

Lemma 6. For n>0,j>2, and < c < 19, 



(38) 



|&,(n,c)| < 



[l + cy 
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Proof. The proof will proceed by induction on both j and n. First, calculations 
show that 



6,(0,c) 



1 ifj = 0, 
if j > 1. 



due to the simple form of Kq{u). Using this calculation together with the recursive 
nature of the coefficients, we also calculate that 

'log(c+l) ifj-0, 



if J > 1, 



and. 



&,(2,c) 



i^2(2 + c) 
log(c+l) 

c + 2 

jic + 2y 



logic 



1=1 



I Vc+l 



ifi = 0, 
if J = l, 

ifj >2. 



The bound claimed in the lemma is therefore clear for n = and n = 1. For the 
case when n = 2, we will need to show that 



c + 2 



c+1 



1=1 ^ 

For one side of the inequality, we have that for < c < 19, 

i-i 

log(c+l)-^ 



1=1 



I Vc+l 



< log(c + 1) - <2<j' * 



c+1 



c+1 



For the other side of the inequality, we must show that 



E 

1=1 



I Vc+l 



- log(c + 1) < J 



c + 2 
c+1 



but here, Grupp and Richert 6, first formula below (4.6)] obtain the superior bound 

'c + 2^^' 



1 / c + 2 



1=1 



I Vc+l 



- log(c + 1) < 



c+1 



Let us therefore assume that n > 3 from now on. Before we can induct on both j 
and n, we need to prove that the bound in (38) holds for j = 2. Here, Grupp and 
Richert W, formula (2.9)] supply us with the useful inequality 

log"(u-n+l) 



(39) 



< K,,{u) < 



This bound clearly holds for n = 0. By induction, when n > 1, we have 



KJu) < 



1 



< 



(n-iy. 
1 



log"-i(t-n + l) 



dt 



log""^(t - n + l)dlog(i - n + 1) 



log"(u - n + 1) 
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since t — n + 1 < t ior n > 1. Thus, since log(c + 1) < 3 for < c < 19, it follows 
from Lemma [sj ([36]) , and ( 39 ) that 



|^'2(n-,c)| 



1 



< 



< 



2(n + c) 
1 

2{n + c) 
1 

2(n + c) 
1 

< — -r max 

- 2(l + c)2 

1 

< — -r max 

- 2(l + c)2 



[n — 1 + c) (n + c) 



{n-l + c) 
log"-2(c+l) 



[n + c) 
log"-i(c+l) 



(n - 2)!(n - 1 + c) ' (n - l)!(n + c) 
log"-^(c+l) log"-i(c+l) 
(n-2)! ' (n-1)! 



< 



(n-2)!' (n-1)! J " (l + c)2' 
To complete the induction, we observe that if j > 3 and n > 3, 

1 



\hj{n,c)\ 



< 



j{n + c) 
1 



{bj-i{n - l,c) - (j - l)6j_i(n,c)} 

4 , , 4 



J 



1/4 (■ ] \ 



< 



(1 + c)^- 



□ 



If one requires a bound for bj(n, c) in a larger range of c values, say 2 < c < C, 
one could probably replace the constant 4 in the lemma above with log(C + 1). A 
bound for c < 19 is more than enough for our purposes. Grupp and Richert [6j 
remarked that the bound in (37) could be extended to hold for < c < 9, but with 
considerably more work. 



6. The Chain of Circles 

In the last section, many facts concerning the power series representations of 
Kn{u) were assembled. This information will be especially useful when combined 
with an idea of Grupp and Richert, known as the Chain of Circles, or Kreisketten- 
verfahren. The method is essentially analytic continuation. To begin, one defines 
the sequence 

and forms the corresponding sequence of power series 



- 1, 



Kn{u\v) = ^6j(n,c^) (u- {n + Cy)y . 
This sequence of power series has the feature that it can be generated recursively. 



The power series for Kn{u\ v) is obtained from Kn{u; j/ — 1) since, using (35), 
(40) bo{n,c„) ^ Kn{n + c^;iy - I), 
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and the rest of the coefficients can be computed using Lenima|4]or Lemmajsj Thus, 
we have a chain of power series representations for Kn{z) that can be continued 
throughout the half plane SRz > n — 1, as seen in Figure below. 




Figure 1. Chain of Circles 

Although the power series Kn{u; v) is valid inside a larger interval, we will restrict 
the representation to the interval 



{u : n + Cu < u < n ^ c^+i} 



to speed the convergence of the series. The sequence c^, as Grupp and Richert 
point out, strikes a balance between the number of power series needed to cover a 
fixed u value, and the convergence rate of each of those power series. Finally, we 
have obtained a useful decomposition of Kn{u), given by 



(41) 



Kniu) 



OO 

E 

1^=0 



where Xi/(u) is the characteristic function of the interval I^. Now, for numerical 
purposes, we will truncate each of these power series to, say, N . Actually, for our 
purposes we will eventually take N = 80. In the first circle, Kn{u]Q) will suffer 
only from the truncation. However, in the next circle, Kn{u; 1) will not only be 
truncated, but the coefficients will be approximates of the actual coefficients due 
to the recursive nature of bo{n,ci) — Kn{n + ci;0). Controlling the error that 
propagates will therefore require some work. To make our discussion more precise, 
let us define 



(42) 



Kn{u; v) 



N 

E 



bj{n,c,y) {u- {n + c„)y 



The coefficients bj{n, c^) will be generated in exactly the same fashion as 6j(n, c,^) 
using (40l and Lcmmajij When ly = we have bj{n,co) — bj{n,co), for j < N. 
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However, the bj{ n, c^) will be approximates of the actual coefficients bj{n,c^) for 
v > 1 due to ( 40 1 . More specifically, we define 

and 



1 ifj=0, 
ifj>0. 



bj{n, Ci,) = < 



Kn{n + c^) 
(-1) 



if j = 0. 



j{n + c„y 



1=0 







if j > N. 



The following lemma of Grupp and Richert tells us that the error between the 
coefficients bj{n,c^) and bj(n,Ci,), for 1 < j < N , can be obtained from the corre- 
sponding error when j = 0. 



Lemma 7. //, for a fixed c> 0, we have 

bo{n,c) - bo{n,c) 

then, for 0<j<N, 



bj{n, c) — bj{n, c) 



< 



s 



(2 + !)■'"■ 



Proof. This is proved by induction on n. It is vacuously true for rt = and n = 1 
since in those cases we will take bj{n, c) = bj{n, c). Now, by induction. 



j-i 



~ 1 I ~ 

bj{n,c) -bj{n,c) < j^^_^^y X! 1'^) ~ bi{n - l,c) (n + c)' 



j-i 



j{n + cy + 



1 



E 



3-1'' 



since the terms in this last sum arc all bounded above by one. 
Following Grupp and Richert, we prove 

Lemma 8. IfO<j<N,v>l, and Ci, < 19, then 

1 M,_i 
{2+fy 2- ' 



□ 



(43) 

where 
(44) 



bj{n, c^) - bj{n, c^) 



< 



7 + ci 



4n 2 



1 /3 



i=0 



3 V2 
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Proof. We are going to use Lemma [7] to establish that for i' > 1, 



(45) 



< 



The proof will proceed by induction on i/. For = 0, we will take bj{n,co) 
bj{n,CQ). Therefore, when v = 1, we use Lemma[6]and (40) to observe that 

bo{n,ci)-bo{n,ci)^ < ^ \bj{n,co)\ {ci - cqY < ^ < ^ ("^T^ 

j>N ^ 

Hence, by induction, the difference 

bo{n,c^) - bo{n,c^) = Kn{n + Cy-.v ~ I) ~ Kn{n + c^;v - I) 
is at most 

\bj{n,Ci,^i) -bj{n,Cy^i) {cy - Cy^iY + ^ \bj{n,Cu^i)\ {c„ - c^-iY . 

j<N j>N 

Using the inductive hypothesis together with Lemma [6j this series is bounded by 



1 + E 



OO / \ J 



3=0 



2+ V 



2N 



v-2 

n 

/=0 



7 + Q 



n 

(=0 



7 + Q 



The induction is complete, and from Lemma [Tj 



bj{n,c^) - bj{n,Cy) 



< 



1 / 3 



(2 + ^r2-ilU 



□ 



Now that we have good control of the coefficients bj{n, c^), we prove the following 
bound concerning the error between Kn{u; v) and Kn{u; v). 

Lemma 9. If n > 0, ly > 0, < 19, and N > 2, then 



(46) 



< 



2N ' 



where is as in (44). 



Proof. The proof will proceed by induction on v. When z/ = 0, we will take 
bj{n,Co) = bj{n,cn), so 

Kn{u;0)~Kn{u;0) = ^ 6j(n, co)(u - (n + cq))^'. 

j>N 

Thus, using (|38|, 



2N - 2N 



~ I x-^ 4 Mn 

Kr,{u-{)) ^ Kr,{u-Q)\ < \b,{n,co)\{ci-coy < ^ < 

]>N 



For V > 1 , we use ( 38 1 and Lemma jsj since 
^ \bj{n,c^) -bj{n,Cy)^{c^+i 

j<N 



Kn{u; v) ~ Kn{u; v) 



is at most 



^ \bj{n,Cy) \ [cy+i - c, 



V } 1 



j>N 
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which is bounded by 



E 



M^_i ^ / c^+i ~cA\ 4 / ^ / c^+i - V \ Af^_i M, 



< 



2^ ^ \ 2 + ^ J 2^-1 2+% ; 2^ 2^ 



7. Generalizing to Integral A 7^ 



□ 



When considering integral A 7^ 0, one is faced with the problem of understanding 
the coefficients of the power series representation 

00 

(47) Kn{u, A) = ^ b,{n, c, X){u - {n + c)y , 

3=0 



again valid inside | u — (n + c) | < 1 +c, by the same reasoning as in ( 34 1 . The critical 
observation here is that repeated applications of Lemma [3] can be used to write the 
bj{n, c, A) in terms of bj{n, c, 0) = bj{n, c). Thus, to generate these coefficients, one 
can use the fact that 

(48) bo{n,c,X) = K„{n + c,X), 
and, for j ^ 0, 

(49) bj{n,c,\) = -bj^i{n,c,X~l). 

The analytic continuation technique of Grupp and Richert will be carried out similar 
to the case when A = 0. As before, these power series will be chained together to 
generate expansions throughout the interval u > n. Thus, one defines 

CO 

Kniu,X;v) =^ ^bj{n,c^,X) (u - {n + c^)y , 
j=o 

each one valid inside the interval Ii, = {u : n + c,y < u < n + c^+i}. This sequence 
of power series can be generated recursively. The power series for Kn{u, X;^) is 
obtained from if„(u. A; z/ — 1) since 

(50) bo{n,c^,X) ^ Knin + c„,X;iy - 1). 

This is precisely how the power series expansions are chained together. The prob- 
lem, of course, is that we will have to settle for an approximation to Kn{n+c,^, A; v— 
1), as this value will be obtained by a truncated power series expansion. The series 
are related to the Kn{u, X) functions via the decomposition. 



(51) Kn{u,X)^Y.X.{u)Kn{u,X:v) 



00 

u=0 



where x,y(M) is the characteristic function of the interval /jy. Of course, we make the 
definition Kn(u, 0; v) — if„(u; v). We produce power series that represent Kn{u, A) 
in various intervals. We will truncate these series for numerical purposes, and hence 
define 

N 

(52) Kn{u,X;v) = '^bj{n,c^,X) (u - {n + c^)y . 
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The coefficients 6j(n, c^, A) are defined by 



if < j < A, 
5,(0,co,A) = <( 1 ifj-A, 
ifj>A, 



and 



Kn{n + Ci,,X- I) if j = 0, 

hj{n,Cu,X) = { c,,, A - 1) if < j < iV, 

if j > N. 

When V = Q we have that co,A) = 6j(n, co,A), for j < N. As before, the 
bj(n, Ci,, X) will be approximates of bj{n, c^, X) for i' > 1 due to (50). In any case, 
we proceed as in (51 1 and set 

oo 

Kn{u, A) = ^ x>y{u)Kn{u, A; ly). 



The purpose of this section is to bound the error between Kn{u, A) and Kn{u, A). 
Thus, we prove 

Lemma 10. // n > 0, i/ > 0, < A < iV, c^, < 19, and N >2, then 



(53) 

where q = and 



Kn{u,X;v) - Kn{u,X;v) 



< 



XlM,,x 

2N-X ' 



(54) 



1/ 1 /3\'' 



k=0 



fe=0 



Proof. The proof will proceed by induction on both v and A. The case A = has 
already been shown in Lenima[9| When = 0, we will take bj {n, cq, A) = bj (n, cq. A), 
so if < A < iV, we can make repeated use of (49) to see that 

Kn{u,X;0) - Kn{u,X;0) = ^ bj{n, cq, X){u - {n + co)y 

j>N 

A A- 1 1 



E 

j>N 



j i - 1 j 



j-^bj^x{n,ca,0){u - (n + co))^ 



Thus, using p§ 



Kniu, A; 0) - Kn{u, A; 0) < A!(ci - co)^ ^ \bj-x{n, co, 0)| (ci - co)^-^ 



< A!(ci-co)' 



2JV-A 



< 



A!(ci-co)H/o,o A!Afo,A 



2N-X 



2N-X 



We have shown that (531 holds for A = 0. To prove (53), observe that Kn{u, A; v) — 
Kn{u, A; v) can be rewritten using Lemma [sj and (52) as 

Knin + C;,,X;iy-l)-Kn{n + c^,X;iy-l) + X / Kn{t, X-l;iy) - Kn{t, X~l;iy)dt, 
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where K is K with N replaced by A'^ — 1. The first two terms above correspond 



to the j = term of the power series expansion. Finally, the bound in ( 53 ) follows 
since 



OJV-A ' ' / OAT-A 

□ 



Although the presence of the A! term in ( 53 ) looks menacing, we plan on taking 
A < 10. In addition, we will take N to be much larger than A, say N = 80, so the 
error will still be well under control. In the next section, we will apply this theorem 
to approximate j'^. 

8. Approximating j'i^{u) in the Main Computation 
Recall the integrals 

(55) Ii = / P(wf j'^{u~w)dw, 

{P{w)-P{w-t)f jj',{u-w)dw, 
(57) ^3=/ / {P{w)-P{w-t)f -j',{u-w)dw, 



JO 

1 

2 



t 



1 '•1 

t 



(58) T4= / / P{wY -j'^{u-w)dw. 



If Ii — ^ — KI4 > 0, then a positive lower bound for S^A,!^, z) is obtained. 
To compute these integrals, define 

(59) Xi = / P{wfJ'^{u-w)dw, 

Jo 

(P{w)-P{w-t)fjf^{u-w)dw, 

(61) 13= / / {P{w)-P{w~t)f -j'^iu-w)dw, 

Ja Ja 

~ dt ~ 

(62) Z4=/ / P{w f -fAu-w)dw, 



t 



where 



^ ' 0<n<u 

Set / = Ii — ^ k2!^3 — KI4, and / = Zi — ^ — KI4. Naturally, the integral 
/ approximates /. The computations below are performed with = 80 in the 
definition of K„{u, n — 1). The error between / and / is bounded using Lemma 
in the last column. 



10 
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Table 2. Sifting Limit Calculations 



K 




u 


a 


1 


Error 


Z 


4.5iD 


1. toiii 




Z.9 X iO 


a o \v 1A — 23 
6.0 X iO 


Q 
O 


D.ozU 


I. /OUi 


U.ZDz /Oi 


K /I 1 A — 6 
0.4 X iU 


0.0 X iU 


4 


8.522 


3.7611 


0.260302 


2.3 X 10-^ 


1.2 X 10-20 


5 


10.523 


4.7617 


0.258785 


4.5 X 10"^ 


2.3 X 10-19 


6 


12.524 


5.7621 


0.257739 


6.7 X 10-5 


4.9 X 10-1^ 


7 


14.524 


6.7623 


0.256929 


2.2 X 10-5 


1.2 X 10-1^ 


8 


16.524 


7.76247 


0.256318 


9.3 X 10-^ 


3.9 X 10-15 


9 


18.525 


8.7627 


0.255870 


6.5 X 10-5 


1.5 X 10-13 


10 


20.525 


9.7628 


0.255468 


4.8 X 10-5 


6.7 X 10-12 



These computations verify the values appearing in Tablenlfor /J^ given by the A^A" 
sieve. These calculations conclude the proof of Theorern^ 

I would like to thank my advisor, Sid Graham, for supervising this work and 
taking the time to carefully read through this manuscript. 
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